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Abstract 



r^ ' The heat operator with a general multisohton potential is considered and its ex- 

tended resolvent, depending on a parameter g £ R^ is derived. Its boundedness prop- 
erties in all variables and its discontinuities in the parameter q are given. As the result, 

K^ the Green's functions are introduced and their properties are studied in detail. 
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The Kadomtsev-Petviashvili equation in its version called KPII 

(U4 -6MMa;i +Wxia;ixi)xi = -3Mx2:C2) (1-1) | KPII 



where u = u(x,t), x = (xi,X2) and subscripts xi, X2 and t denote partial derivatives, is 
a (2-|-l)-diniensional generalization of the celebrated Korteweg-de Vries (KdV) equation. 
k^ ' The KPII equations, originally derived as a model for small-amplitude, long- wavelength, 

^ . weakly two-dimensional waves in a weakly dispersive medium [1], have been known to be 

integrable since the beginning of the 1970s [2,3], and can be considered as a prototypical 
(2-f l)-dimensional integrable equation. 

The KPII equation is integrable via its association with the operator 

C{x, dx) = -9x2 + dl^ - u{x), (1.2) heatop 

which define the well known equation of heat conduction, or heat equation for short. The 
spectral theory of the operator (1.2) was developed in [4-7] in the case of a real potential 
u{x) rapidly decaying at spatial infinity, which, however, is not the most interesting case, 
since the KPII equation was just proposed in [1] in order to deal with two dimensional weak 
transverse perturbation of the one soliton solution of the KdV. The main difficulty in the 
study of this problem is due to the fact that the soliton solutions of the KPII equations 
are not decaying at space infinity and have a ray behavior on the a;-plane (see, e.g., [8-11]). 
Correspondingly, the integral equations defining the Jost solutions in this case are senseless, 
as their kernels do not exist. A spectral theory of the KPII equation that also includes 
solitons has to be build, as already successfully done for the KPI equation [12]. Following a 



framework of this article we develop a generalization of the standard 1ST so called "scattering 
on nontrivial background," i.e., we consider a potential 



u{x) = u{x) + u' (x) , X — {xi,X2), 



(1.3) 



where u{x) is some multisoliton potential and u'(x) is a smooth, rapidly enough decaying 
function of its variables, that can be considered a perturbation of the soliton potential. 
Correspondingly, the Jost solution of the operator C with potential u can be defined as the 
solution of the following integral equation 



$(a;,k) = $(a;,k) + / dy ^(a;, y, k)u'(y)$(y, k). 



where k g C is a spectral parameter, $(a;, k) is the Jost solution of the operator C with 
multisoliton potential u and Q{x,y,'k.) is its total Green's function, i.e.. 



(-9., + dl - u{x)) g{x, x\ k) = S{x - x'). 



(1.5) 



In order to be useful for developing the 1ST on the solitonic background, the Green's function 
must obey condition of boundedness, i.e., function 



G{x, x', k) = e'k(:ri-^'i)+k'(^2-4) g(^x, x' , k) 



(1.6) 



must bounded with respect to the variables x,x' G M and k G C and has finite limits 
at infinity. In this framework it was possible to develop the inverse scattering transform 
for a solution describing a perturbation of the one soliton solution [13]. However, the case 
of any number of solitons is still open. For this case, in [14] we derived a total Green's 
function, which is a natural generalization of the Green's function for the case of a decaying 
potential. But as it was shown in [13] and [12], in order to control the singularities of the 
Jost solutions we need some other Green's functions. Thus, following [12], in order to deal 
with a heat operator with a generic multisoliton potential (for some multisoliton solutions 
this was done in [15]), we introduce its extended resolvent, which is an object more generic 
than the Green's function. Exactly, we introduce a two dimensional real spectral parameter 
q = (91,(72) and consider the extended Lax operator 



C{x, d^+q) = -d^^ - (72 + {dxi + qif - u{x), 



and, then, the extended resolvent of the heat operator (1.2) is defined as the tempered 
distribution M{x, x'; q) with respect to all its six variables x, x', q that satisfies the differential 
equations 

C{x, dx + q)Mx, x'; q) = C^ix', d^' + q)M{x, x'; g) = <5(x - x'), (1.8) 

where C'^{x, dx) is the dual oi C{x,dx)- 

The extended resolvent can be considered the generating functional of the different 
Green's functions of the /^-operator (1.2). Indeed, let us introduce 



M(x,x'; q) = e 



— pii^-x') 



M(x,x';q), 



(1.9) 



that, of course, is not necessarily a tempered distribution. Nevertheless, it is easy to see 
that 

C{x,dx)M{x,x';q) ^ C^i^x ,dx')M{x,x'; q) = 6{x-x'). (1.10) 



(1.4) tphi 



green 



Green 



(1.7) heatop2 



risolvl 



Mhat 



MhatGreen 



In particular, the total Green's function is given by means of the reduction 



g{x,x',k) ^ M{x,x';q) 



9i=ka 
92=k| - 



(1.11) 



The construction of the distribution M{x,x';q) is the subject of this article. In [16] it 
was shown that the heat operator (1.2) with a multisoliton potential u{x) can have left or 
right annihilators in polygonal regions of the g-plane, where, consequently, the resolvent 
cannot exists. Here we give an explicit expression for the extended resolvent and prove that 
outside some special polygon region (the same as was given in [16]) it exists as a tempered 
distribution and obeys (1.8). We demonstrate that in this region reduction (1.11) is always 
possible for an arbitrary k e C and it gives the total Green's function derived in [14]. Thus 
the condition that M{x,x';q) is a distribution leads to the boundedness property (1.6). 
According to the procedure followed in the case of the nonstationary Schrodinger operator 
in [12], the description of the discontinuities of the total Green's function on the complex 
plane of the spectral parameter requires some special reductions of the extended resolvent, 
i.e., auxiliary Green's functions. Thus we conclude this article by a detailed study of these 
discontinuities. 



GM 



2 Heat operator with multisoliton potential and its Jost 
solutions 

Soliton potentials (see [15] and [8-11,16,17] for details) are labeled by two numbers (topo- 
logical charges) Na and Nb, which obey condition 



iV„,7Vh>l. 



Let 



so that N >2. We introduce the N real parameters 

Kl < K2 < . . . < KU, 
and the functions 



Let 



K^{x) ^ n^x-^+K^^x^, n = l,...,A/'. 
e^(-)=diag{e^"(-)}^=, 



be a diagonal N xM matrix, let V he a. N xNt, constant matrix and V be an "incomplete 
Vandermonde matrix," i.e., the Nf, x Af matrix 





(I ... 1 \ 


v = 


Ki ... K_A/ 




\ 1^1 ■ ■ ■ I^J\f / 


Then, the soliton potential is given by 


V 


(x) = -2dllogT{x), 



(2.1) nanb 




(2.2) Nnanb 




(2.3) kappas 


(2.4) Kn 

(2.5) eK 



(2.6) 



(2.7) ^ 



where the T-function can be expressed as 

t{x) ==dct(Ve^(=")l?). 



(2.8) 



tau 



For the Jost and dual Jost solutions (solutions, respectively, of the heat operator (1.2) and 
its dual) we have 



$(x,k) =e-*''="i-'^'^^x(a;,k), 
*(x,k) = e*''^i+'''^^e(2;,k), 



(2.9) 
(2.10) 



Phi 



Psi 



where 



X(x,k) = 



C(x,k) 



$(^; k) 

t{x) 

T^l,{x,\i.) 

t{x) 



(2.11) 


symPhi 






(2.12) 


symPsi 



with (Miwa shift) 

. (x, k) = det (V e^("=) (k + i k) V) , t* (x, k) = det [ V - 



K{x) 



T*i 



+ zk 



V 



(2.13) 



tauk: 1 



and 



K + i k = diagJKn + i k}„^;^ 



(2.14) kappadiag 



In order to study the properties of the potential and the Jost solutions, it is convenient 
to use the representations for the r-functions, which follow by the Binet-Cauchy formula 
for the determinant of a product of matrices, i.e.. 



^f 



Nt 



^(^) = 1^ E T^i{nr})V{{n.})Y[e^-'^^\ 



ml 



{«,} = ! 



1=1 



H 



^^^'^^^TV^ E V{{m.])V{{m.})X{{nm,+i\^)e^-^^^\ 



{mi} = l 

AT 



^(->k) = -rr|^ E ^(KDnKDlI 



1=1 



Nb\T{x) 



{",} = ! 



1 = 1 



^ni n" ^ K 



where we used notations 



(2.15) 
(2.16) 
(2.17) 



taufl 



Phikl 



Psikl 



F({n,}) = det 



.N,,~-l 



T^ni,l 



V{{n^}) = det 



.Nt-l 

V. 



l<i<j<Nb 



ni,Nt 



v„ 



for the maximal minors of matrices V and 2? and where 

{rrii} = {mi,...,mjvj, {nj = {ni,...,nArJ 



(2.18) 



(2.19) 



Dol 



(2.20) |1^ 



stand for non ordered sets of Nt, indices from the interval 1, . . . ,Af. We recall that the 
maximal minors of a matrix satisfy the Pliicker relation, i.e., for any subsets {rrii} and {rii} 
of indices running from 1 to J\f and arbitrary j G {1, . . . , Nf,} 



V{{m,})V{{n,}) 



Nb 



y^ T>{mi, . . . , TOs_i, rij , TOs+i, . . . , ttiatJ P(ni, . . . , rij^i^ms, rij+i, 



,nN, 



(2.21) 



plu 



Notice that the only x-dependent terms in (2.15), (2.16), and (2.17) are exponents of 
sums of linear functions (2.4). Correspondingly, the asymptotic behavior of the function 
t{x) and of the potential has a sectorial structure on the x-planc. In order to specify these 
sectors at x — > cxd we introduce the ray directions 



xi + (ku + 'in+Nt)x2 bounded 



n = 1,...,A/', 



(2.22) [TF] 



where we assume that the indices are defined mod A/", so that thanks to (2.2), say, n + Ni, — 
n — Na for n > Na- Thus there are Na rays in the direction a;2 — > — oo and Ni, rays in the 
direction X2 — s- +oo. The sector an is swept out by rotating anticlockwise the ray r„ up 
to the ray r„+i. These sectors are nonintersecting and cover the whole x-plane with the 
exception of rays. In [9] we proved that the the leading exponents of t(x) when x — > cxd 
are the exponents exp(^"J^ '" Ki{x)) , each being the leading one in the corresponding cr„ 
sector of the x-plane. More exactly, if the coefficients 



z„ = V{n,...,n + Nb - l)V{n, . . . ,n + Nb - 1), 



(2.23) [^ 



are different from zero for all n = 1, . . . , A/" (again with indices defined mod A/") the function 
t{x) has, along rays and inside sectors, the following asymptotic behavior 



-{x) - (z„ + z„+ie^"^+"(-)-^"(-) + o(l)) 



-(x) := (z„ + o(l)) cxp ^ Ki{x) 

\ l=n J 



fn+N^,-\ 



exp 



^ K,{x)Y (2.24) 



(2.25) 



Regularity of the potential u{x) on the x-plane is equivalent to the absence of zeroes of 
t(x). It is clear that it is enough to impose the condition that the matrix T) is Totally Non 
Negative (TNN), i.e., that 



P(ni 



, '^AfJ > 0, for aU 1 < ni < . . . < rijVi, < M . 



However, sufficient conditions on the matrix 2? for the regularity of the potential are un- 
known. On the other side, from (2.24) and (2.25) it follows directly that it is sufficient to 
require that 

z„ > (2.27) 

for having nonsingular asymptotics of the potential. 

We also mention that the functions x(x,k) and ^(x,k) have bounded asymptotics on 
the x-planes because the x-dependent exponents enter in denominators and numerators of 
expressions (2.16) and (2.17) with coefficients proportional to 2?({rii}). This means that the 
leading asymptotic behavior of the denominators of the functions x(a;, k) and ^(x, k) on the 



3:232 



tauasympt 



(2.26) 



znO 



x-plane is not weaker than the behavior of their numerators. For more details see [8,9,16,17], 
where the same notations have been used. 

We need in the following also the values xi^j *'*«) of x(x, k) at k = iK„ and the residues 
C„(a;) of C(a;,k) at k = iK„. From (2.9), (2.10) and (2.16), (2.17) we have 



X(x, iK„ 



(-1) 



N, 



Nb\T{x) 

1 



iNb\T{x) 



J2 V{{m.})V{{m,},n)Yle''-'^^\ (2.28) 

{mi} = l 1=1 

J2 P({n.})^<5„^„(-l)^~V(ni,...,^,...,n^jn^'^"'''^ (2-29) 



where {{mi},n} ~ {mi, . . . ,mN^,n}, hat over rij denotes that this index is omitted and 
where the (5„.„ Kronecker symbol in the r.h.s. of the last formula is due to the fact that the 
residues of the terms in the sum are nonzero only when some rij — n. 

Taking into account the analyticity properties of x(a;, k) and ^(x,k) in (2.16), (2.17) 
their product can be written in terms of the values x(a;,iK„) and ^n{x) as follows 



A/" 



X 



{x,k)ax',k) = l + J2 



X{x,iKn)^„{x') 



n=\ 






(2.30) 



which also will be useful in the following. In [16] we demonstrated that the Jost solutions 
obey the Hirota bilinear identity 



M 



^$(x,iK„)*„(x') ==0, 



(2.31) 



n=\ 



where in analogy to (2.29) ^„(x) denotes the residue of *(x,k) 



chixi 



sumPhiPsi 



3 Extended resolvent M(x, x'\ q) 

Here we prove that the resolvent M{x,x';q), i.e., a tempered distribution of its 6 real 
variables that obeys (1-8), can be written as a sum of a continuous (in some sense) and a 
discrete part 

M(x,x';g) =Mc(x,x';g)+Md(x,x';q), (3.1) \k\ 

whose definitions and properties are given below. We also specify the infinite region on 
the (/-plane where the property that M(x,x';q) is a tempered distribution holds, and we 
prove, in fact, that M{x,x';q) is bounded for q in this region. For both terms we use the 
hat-kernels notation introduced in (1.9) 



Mc(x, x'; q) =. e-«("-"')Me(x, x'; q), M^x, x';q)^ e-'^'^^-^">M4x, x'; q). (3.2) 

So, first, we define 

sgn(x2 - 4) 



Mc(x,x';g) 



27r 



da e{{q2 +a^ - ql){x2 - 4)) ^ 



X $(a;, a + iqi)'^{x , a + iqi), 



(3.3) 



where the Jost and dual Jost solutions $(a;,k) and '^{x^'k) are defined in (2.9)-(2.12), j da 
denotes integration along the real axis, and 6 is the step function of its argument. Properties 
of Mc are given by the following lemma. 



McMDelta 



20 



Lemma 3.1 The integral in the r.h.s. of (3.3) exists and the function Mc given in (3.2) is 
a bounded function of its arguments for all x, a;', g G K. and has finite limits at infinity. 

Proof. Thanks to (2.9), (2.10), and (3.2) we can write that 

M^{x,x' -q) ^ Mo{x,x' -q)- 

,SgTL{x2 - x'2) f, 9{{q2 + a'^ -qf){x2- X2)) 



y2x(x,iKn)£,n{x') 1 -da- 



■X 



X g-ia(xi-x[+2qi(x2-x'2))-{q2+a^-qi)(x2-x'2) j-g ^\ 

where we used (2.30) for the product of x(x, k)^(a;', k) and where Mo{x, x'; q) is the extended 
resolvent of the operator (1.7) in the case of the zero potential 

Ma{x,x'- q) ^ _ ^Sn{x2^x'2) J ^^ ^^^^^ ^ ^2 _ ^2)(^^ _ ^^, )) X 

X g-*0!(2:i-2:i+2(iri(a;2-X2))-((j2+a^-g?)(2;2-a;2) /g r^\ 

Then the statement of the lemma follows directly thanks to the exponentially decreasing 
factors in (3.4) and (3.5). ■ 

Applying the heat operator (1.2) to M(x,x' ,q)c in (3.3) we get 

C{x, d^)M^{x, x'- q) = ^""^^^ '^^' fds $(x, s + %)*(a;', s + iqx). 

The integral in the r.h.s. can be explicitly computed thanks to (2.9) and (2.10) and after 
inserting (2.30) in it, we get 

5{x-2 — x' ) 
^ ^ lds^{x,s + iqi)'^{x',s + iqi) = 5{x - x') - 5{x2- X2)p{x,x',qi)+ 



20:1 



20:2 



g32 



Phat 



2tt 

+ iS{x2 - x'2)e{x[ - xi) Y^ $(x, iK„)*„(x'), (3.6) 

n=l 

where the last term annihilates thanks to (2.31) and where we denoted 

AT 

p{x, x', gi) = z ^ 9{qi - Hn)^{x, iK„)*„(x')- (3.7) 

n=l 

Notice that, p(x,x' ,qi) does not belong to the space of Schwartz distributions, as it can 
have exponential growth in some directions on the x-plane, [16]. Thanks to (2.31) it can 
also be rewritten in any of the following forms 

. A/" A/" 

p(x, x', 9l) = -^ X! ^g^('?i ^ '^n)^{x, iKn)^n{x') = -« ^ ^('«« - qi)^{x, ««„)*„ (x'), (3.8) symPhat 

n— 1 n— 1 

and 

p{x,x' ,qi) — 0, for all qi^[Ki,KAf]. (3-9) | 5.3 | 

Notice also that, just by definition, 

C{x,d^)p{x,x',qi) = for any g e M^ (3.10) |5.31 



Thus we define the second term in (3.1) as (see also (3.2) 

Md{x,x';q) = tO{±{x2 - X2)) p{x , x' , qi) , 

so that for any choice of sign in the r.h.s. we have thanks to (3.10) 
C{x,da;)Md{x,x';q) = S{x2 - x'2)p{x,x'qi), 

that together with (3.6) proves that 

M{x, x'; q) = Mc(x, x'; q) + M^ix, x'; q), 



(3.11) 



(3.12) 



obeys (1.10), or M{x,x';q) given in (3.1) obeys (1.8) thanks to (1.9). Proof of the second 
equalities in (1.8) and (1.10) is analogous. Let us notice that thanks to (3.9) 



Md{x,x';q) = 0, for aU qi ^ [ki,k^]. 



(3.13) 



Thus in order to prove that M{x, x'; q) is the extended resolvent, we need to prove that 
Md(x, x'; q) belongs to the class of tempered distributions and to specify the choice of signs in 
(3.11). For this aim we need to make explicit the dependence of Md(x, x'; q) on its variables. 
First, we consider function p{x,x',qi) defined in (3.7). Inserting in the r.h.s. (2.28) and 
(2.29), thanks to antisymmetry of minors of matrices T) and V (see (2.18) and (2.19)) and 
after summing over n, we get 



p{x,x',qi) 



-1 



A/" 



M 



IS Y. Y. ^^{{m^})v{{n,})e{q^- K y. 



NbUNh - 1)It(x)t(x') 

^ ^ ^ ^ ^ ^ {m,} = l {«;} = ! 

X V{{m,},n^jV{ni, . . . ,njy^_^)x 
xexpiY, ^^i (^) + ^"iv, (^) + Y ^"' (^') 



MDelta 



M+MDelta 



5.32 



(3.14) 



Next, we substitute the r.h.s. of (2.21) with j = Nb for the product of minors of matrix 
T) and exchange mg O nN^, for s ~ 1, . . . , N^. Notice that under this transformation the 
first Vandermonde determinant changes sign, while the second Vandermondc determinant 
is unchanged, as well as the exponent. Thus, we have 



p{x,x',qi) 



X e{q-^ - KmjV{{mi},nNb)V{ni, . . . ,nNb^i)x 
X exp ^ K^, (x) + Kn^ (x) + Y, Kn, {x) . 



(3.15) 



a=i 



i=i 



Exchanging now tti^ o m^ we get, summing over s, Nt equal terms. Finally, we multiply 
(3.14) by Nb, sum up with (3.15) and we divide this sum by A^f, + 1 getting 



p{x,x',qi) = 



1 



{{Nb-iy.)HNf, + iMx)T{x 



r, T. T. 2?({mJ)P(K})) 



{mi} = l {ni} = l 

X [0{qi - K„i„ ) - 9{qi - Knj^ )]V{{Tni},nNjV{ni,n2, ■ ■ ■ ,nNt-i)x 



/ Nt Nt-1 \ 

X exp ^ X„, (x) + Kn^^ (x) + J2 Kn, {x') , (3.16) 

so that for k„ < (ji < k„+i, n = 1, . . . , A/"— 1 

p{x,x',qi) ^ -— ^...2fAT , i\ t \ ^ '\ Yl V{{mi})V{{n,})x 

{{Nb - iy.y{Nb + 1}t{x)t{x') ra,^n~l^...,M 



5.1 



\l^m ?,T i'^ 



Nt'"'"Nh 



]D[k„,K„+i] 

X [e{qi - Km„J - 0{qi - '«n„J]^({mi},njvjT^(ni,n2,...,nAr^_i)x 



Af, Af,-1 



exp Y. ^™< (^) + ^««. (^) + H ^«< (^') ' (3-17) 



5.2 



i=l 1=1 



where the summation runs on all mi, . . . , rriM,, ni, . . . , um,, from 1 to A/" such that the interval 
[Kn, Kn+i] belongs to the interval [«;„„ , k„„ ]. Notice that thanks to V{{mi}, un^) this sum 
does not contains terms with rrijVi, ~ fijVi, • 

Thanks to (3.2) and (3.11) we have that for Kn < qi < Kn+i, n — I,. . . ,7V — 1 

^<^(^'^ ''^) ^ ((I'-W^7^^, + ;MIMlo ^^^„g^^^^^ 2.(WJ)P({nJ)x 

X [0{qi - Km„J - 0((?l - K„„J]F({TOi},njvjT^(ni,n2,...,nArj,_i)x 
/ Wt A^b-l \ 

x exp ^ K^, (x) + Kn^^ (x) + J2 Kn, {x') (3.18) \&\ 

\ 1=1 1=1 J 

with the same condition on summation. This representation for Md(x,x';g) gives another 
proof of (3.13). 

In [16] we demonstrated that the extended operator (1.7) can have annihilators when q 
belongs to some polygons on the g-plane. Thus, the inverse of this operator (the extended 
resolvent) cannot exist for any value of q. We introduce on the q-plane the polygon V 
inscripted in the parabola q2 = q\ of the (/-plane (see Fig. 3), with vertices at the points 
(k„, k^), for n = 1, . . . , A/", whose characteristic function is given by 

A/"-! 

<q) = Y. [^('^1 - '*"+i) - ^('^i - «")][^(9":n+i) - S{qiN)], (3.19) 



19 



where 

qrnn = <2'2 - (^m + Kn)?! + Km^^n- (3.20) 

It is obvious that this polygon divides the strip m < qi < njs/ on the g-plane in two 
disconnected parts. Moreover, this polygon consists of substripes given by subsequent values 
of k's as follows 

K„<gi<K„+i, g„,„+i > 0, n= 1,...,7V-1, gi,A/- < 0- (3-21) 

Now, taking (3.13) into account we can prove the following result. 



qrnn 



sub 



f igur e 




Figure 1: Polygon V in the case TV = 5 



lemma2 Lemma 3.2 Let the {Na, Nh}-soliton potential u{x) be such that its t- function (2.8), (2.15) 
obeys condition that on any subset {ni < . . . < N}y\ G {1, • • • , A/"} the ratio 



n 



eKr.,{x) 



t{x) 



(3.22) [IT] 



is bounded for all x and has finite limits at space infinity. Then Md{x, x' ; q) for all q in the 
strip Ki < qi < Kj\f and outside the polygon V , is a hounded function of all its arguments 
including values at infinities, provided that in (3.11) the upper sign is chosen for q above the 
polygon V and the bottom sign for q below the polygon V . 

Proof. Thanks to (3.13) we consider only q belonging to the strip ki < qi < k^. Let us 
denote Zmn — xi + (nm + '«n)a;2- Then by means of the identity 

0{qi - Krn)-0{qi - Kn) = Sgn{z„in - ^-'mn)'^ 

we rewrite (3.18) in the form 

. , ^ T^(±(a;2-4)) 

Md(x,x';g)- + V V 2 2)} 



((iVb-l)!)2(iVb + l)r(x)T(x') 
^ 2?({TO,})P({n,})sgn(z, 

mi ,ni — 1,...,A/' 



mjV(,njVi, -^mjv^rijvj^ 



F({mi},nArjF(ni,n2, . . . ,nAr(,_i)x 

X [6'((gi - Km„J(2;™ „ - 4»„ „„ )) - Q{{q\ - Knr^Jiz^^ „ - z' ))]x 
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X exp ^ X„, (x) + Kn^^ (x) + J2 Kn, {x') - q{x -x')\. (3.23) \j] 

We decompose this representation as a sum of two terms in correspondence to the two 
terms in the forth hne, we replace q{x — x') in the exponential factors with the identity 
q{x-x') = Km{x) - Kmix') + qmnix2 - x'2) + (qi - K„i){zjnn - z'mn)^ where Km (x) is defined 
in (2.4) and qmn in (3.20), and, finally, we choose m = ttin^, n — tln^ in the first term and 
viceversa in the second one. Thus we get 

Md(x,x';(7) ==M(i)(x,x';q) + M(2)(x,x';g), (3.24) | AP54 

where in each substripe k„ < gi < Kn+i, n = 1, ... ,7V 

.Ai), I ^ _ Tg(±(a^2 - 4)) v^ , , . 

M [X,X ,q)- ((^^_;l)|)2(^^_^1) Z^ Sgn^2„„^n«, ^m«,n«JX 



\l^m^T 1^ 



lVu'"-"JV, 



]D[k„,k„+i] 



xe «™«."«.("^ "^V({mJ,nwjF(ni,n2,...,nw,_i)x 

X e((gi - «;„„ )(z„, „ ~z' „ ))e"*''^"''""J*^"«,"iv,-^™«,"«,^ 

V{{m,}) exp I E;IV' Km, (x) + i^„„^ (x) j 

X ^^ ^x 

r(x) 

P({na) cxp I E;IV' ^«, (x') + Km..^ (x') j 

and 

M<-^hxx'-a)- ^^(^(^2-4)) y- , 

^ ' ' ^^ " ((iVb - l)!)2(iV, + 1) ^ , sgnl2™«,n«, 

mi,ni = l,...,M 

X g-'-""."". ^"^""^ V({mJ, n^jF(ni, n2, . . . , nAr,-i) x 

X e((gi - '^«„J(^™„^„„^ - ^:„„^„„J)e"*''^~''""'^*^-"^"~^"^-"^""^^: 

2?({m,})expfE^'i^™,(^)) 
X Y^ '-X 

P({n.})expfE;i\i^«,(^')j 

r(x') 



(3.25) 



(3.26) 



We see that the dependence of these relations on space variables is due to the exponential 
factors and r-functions. On the other side, since the extended resolvent must belong to the 
space of tempered distributions, we have to show that these exponential factors cannot grow 
either thanks to the ^-functions, or thanks to the r-functions in the denominators. 

Let us consider the x-behavior of these two expressions in detail. First, recalling definition 
(3.20) of qmn, let us notice that, for q belonging to the n-th substripe and g„,„+i < 0, i.e., for 
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q in this substripe below the polygon V, or on its bottom border, all other g^^ ^„„ involved 
in the summation are nonpositive. As well, for qi^j\f > 0, i.e., for q in this stripe above the 
polygon T', or on its upper border, then all other gm„ .„„ involved in the summation are 
nonnegative. Therefore, if one is choosing the signs in the r.h.s. of (3.25) and (3.26) as 
indicated in the Lemma, the exponents in the second lines are decaying or bounded when 
X or g tends to infinity. The exponents in the third lines of the r.h.s. of (3.25) and (3.26) 
are decaying or at least not growing thanks to the (^-functions in the second lines. So we 
have to check the behavior with respect to x and x' of the last two lines of these equations. 
About (3.26) the situation is trivial. The exponents in the forth and fifth lines have the same 
coefficient (minor of V) as in t(x) and t(x') in the denominator and, therefore, these terms 
are bounded when x and x' are growing. Situation with (3.25) is more involved. Let us 
consider the term in the forth line. If the minor 2?({TOi}) in the numerator is different from 
zero and T>{mi, . . . , tojV(,_i, riNb) 7^ 0, then the same exponent as in the numerator is present 
in r(x) in the denominator and the ratio is bounded. However, if ^(toi, . . . , nipf^-i, un^) — 
such exponent is not involved in t{x) and in the direction where it is the leading one (if 
such direction exists) the ratio is growing at large space. The same is valid for the term in 
the fifth line of (3.25). Thus boundedness of these ratios and then of the whole expression 
(3.25) is guaranteed by the condition (3.22) of the lemma. ■ 

Remark 3.1 It is clear that condition (3.22) of this lemma is enough for the validity of 
its statement, hut not necessary, and boundedness of Md{x,x';q) requires additional study. 
Nevertheless, it is already clear that the case of a Totally Positive (TP) matrix D guarantees 
the implementation of condition (3.22). If instead of TP we impose on the matrix T> condi- 
tions (2.27), then all leading exponents are involved in t{x), as we mentioned in discussion 
of (2.24) and (2.25). Thus, (3.22) holds if t{x) has no zeroes in the finite domain. To avoid 
this singularities it is enough to impose additionally that the matrix T> is TNN. 

Remark 3.2 Boundedness of M{x,x';q) with respect to the variable q when g G M \V , on 
the boundaries ofV, and in the limits at q-infinity follows from boundedness of Mc{x, x'; q) in 
(3.4) and M^^'{x,x';q), M^'^'{x,x';q) in (3.25) and (3.26). In the next section we consider 
the behavior of the function M(x, x'; q) with respect to q in detail. 

Remark 3.3 Summarizing, we proved that under conditions of Lemma 3.2 the function 
M{x,x']q) is a bounded function of its arguments and has finite asymptotic behavior. This 
means that this function belongs to the class of tempered distributions, i. e, it is the extended 
resolvent of the heat operator C{x, x'; q) for q outside the polygon V . 

4 Properties of the resolvent and Green's functions 

4.1 Extended resolvent inside and outside parabola q^ = qj 

Local properties of the extended resolvent are easier to study in terms of the hat-kernel (1.9). 
Under the special reduction (1.11) this kernel is nothing but the total Green's function 
t/(x,x',k), where k G C is the spectral parameter. In fact, since under this reduction 
(l2 — q\ = ^ k^ is nonpositive, and qijv (see (3.20)) are less or equal to zero, it gives a 
mapping of the exterior of the parabola ^2 = Qi region of the q plane on the complex plane 
of k (precisely, a one-to-two mapping since the reduction depends on | ks^ | and not on ks){). 
Taking (3.13) into account we see that the part of the strip ki < qi < Kj\f outside of the 
parabola is below the polygon V, so we have to choose the bottom sign in (3.11). Thus, 
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thanks to (3.1), (3.3), and (3.7) wc get 
sgn(a;2 - x'2) 



Q{x,x',k) 



2tt 



ia\X2 — X2 



dad(^{a^ — ks,j)(x2 — x'2))^{x,a + ikcj)^(a;',a + ikcj) + 



Af 



Y^ 9{k<^ -K„)$(a;, m„)*„(x'), 



(4.1) I propl 



that is the Green's function derived in [14]. Its boundedness property (see (1.6)) proved 
there follows now from boundedness of the resolvent at infinity. 

The interior part of the parabola, as was already mentioned in discussion of (3.19), 
is divided by the polygon V in the part above the polygon {qi.j\f > 0) and A/" — 1 lenses 
bounded by the parabola and its chords connecting points q — (k„, k^) and q — (k„_|_i, k^+i), 
n = 1, . . . , A/" — 1 (see Fig. 3). All these lenses are below the polygon V. Taking into account 
that for (72 ^ Qi both Mc{x,x' ;q) and Md{x,x' ;q) arc independent of (72 and in analogy to 
(1.11) we introduce 



g~^{x,x\'kQ-) = M{x,x':q) 
g'{x,x'M)^M{x,x'-q) 



<3i=ka 
i3i=ka 



92 > 9l7 92 > (ki + Kj^)qi - KlKA/-, (4.2) 

92 > 111 92 < (Kn + K„+i)gi - K„Kn+l, (4.3) 

n = 1,...,7V. 



Again, thanks to (3.1), (3.3), (3.7) and (3.11) (once chosen the proper sign according to 
Lemma 3.2) we derive that 

g+(a;,x',k3) =- ^^^ ~ ^^' / da$(x,a + ikQ)*(a;',a + ik3)- 



AT 



GM+ 



GM- 



ie{x2 - x'2) Y^ 6'(k3 -K„)$(a;, iK„)*„(x'), 



(4.4) I prop2 



g (x,x',ka) 



(^2 ^2) 
2tt 



|da^x,a + ^k,M^',a + ^k..H 



Af 



+ i&(x: 



2 - 2:2) X! ^C*'^ -~Kn)^{x, m„)*„(x'). 



(4.5) 



prop3 



Boundedness property for these Green's functions follows from the boundedness of the 
extended resolvent proved in Lemma 3.2. Precisely, thanks to (1.9) and (4.2), (4.3) we have 
that 



Jia{x[ x^)+s{x'^ x^) g-^i^x.x' ,\<i.<:s) is bounded for any real s such that 
s > kcj + max{0, {K^f - k<_i){'kQ -ki)}, 



and 



(4.6) I prop4 



gkQ(x'i xi)+s{x'2 x^) g (2; 2;'^ kg) is bounded for any real s such that 

k|<s<k|+(K„+i-ka)(k9-K„), n=l,...,AA-l. (4.7) 



prop5 
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4.2 Discontinuities of the Green's functions 

In order to describe the discontinuities of the considered reductions of the resolvent, it is 
convenient to use for the resolvent representation (3.4), so that, thanks to (2.9), (2.10), (1.9), 
(3.12) and (3.7), (3.11), we can write 

M{x, X ;q) — Mo{x, x ; q) — 

, i sgn{x2 - x'2) f . e{{q2+a^ -qDix^-x'^)) 



-^<i>(.,^.„)*„(.')p%^/dJ 



X 



a + i{qi - Kn) 

X g(gi-K„-*a)(xi-a;0-[("+^9i)' + «^'](^2-4) ± ie{±ix2 - 4))6'(gi - K„) I , (4.8) [^^ 



where by (1.9) and (3.5) 

Moix,x';q) = -':^^^^^-^Jda0{{q2+a'-ql){x2-x'2))x 

and where the sign in the r.h.s. of (4.8) must be chosen in agreement with Lemma 3.2. It is 
easy to see, that outside the parabola (72 = Qi the kernel M{x, x'; q) is a continuous function 
of q, when q^ (k„,k^), n = l,...,M. 

From (4.8), by using the reduction (1.11), we get for the total Green's function given in 
(4.1) the following representation 



prop? 






k'^=k^ 



g-(«„+»k')(:r,-<)-(K^^+k' )(x2-x;,)__^^^^_^_^)ggj^^j^^^__^^-)l (4^^Q) I propel 



where for p{x.x'; qi) we used expression in (3.8). 

Thanks to the above discussion, it is a continuous function of k G C for all k 7^ i«;„. In 
order to study the behavior of the Green's functions nearby these points, we consider first 
the Green's function G+, as given in (4.4). Thanks to (4.8) we readily get 

g+(x,x',k3) = JA5l_^ f rfk^ g-.k'(.,-.;)-k-(.2-4) - Vci>(x, *«:„)*„ (x')x 

{1 r ^ e-(Kr.+ik')(a;i-a;i)-(K^+k'^)(a;2-2;2) j I 

2^ Mk^ l7~i^ + -sgn(k3-K„) ^. 

(4.11) 

It is easy to see that this Green's function is continuous with respect to the variable kg, 
and, moreover, Q'^{x,x' jk^) is independent of kcj. Indeed, differentiating by kcj, we use 
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propS 



that exponents are analytic with respect to the variable k' = k'jj +i kg and rapidly decay 
when k^ tends to infinity, and also that 5k'(k' — zk„)~^ = 7r(5(k^)(5(k'cj — k„). Thus the 
derivatives of the first term and every term in the parenthesis are equal to zero, so that we 
can write that 



^$(x,^«;„)*„(x') / _^e-*"(^i-^i+2''"(^^-^^»-"'(^=^-="^) L (4.12) 



ra=l 



where integrals are understood in the sense of the principal value and that also can be 
calculated explicitly in terms of the hypergeometric functions. 
On the other side, directly by (4.4) and (4.5) we have that 



. AT 

I 



g (x,x',k3) =g+(a;,a;') + -^sgn(kQ-K„)$(a;,iK„)*„(x'). (4.13) 

^ n=l 

So tj^(x, x', kcj) is discontinuous at all kg — k„, and, precisely, 

t/"(x,x',K„ + 0) - ^"(x,x',fi;„ - 0) == i$(x,iK„)*„(x'), n = l,...,7V. (4.14) 

Thus, considering the difference of the expressions given in (4.10) and (4.11) we derive 
that the total Green's function can be written in the form 

g{x, x', k) = erog(a;, x', k) + g^{x, x', k), (4.15) 

where 

—..-4/".- • " 

-|k»| 
1 .^ r g-(Kn+Jk')(a;i-rri)-(K^+k'^)(i;2-2:2) — 1 

+ ^ II *(^' *'*")*« (^') / '^^'r V^- ' (4.16) 

ZTT ^_^ J K —IKn 

-|k!R 
k'a^kg 

the last integral term being regularized by a subtraction of 1 in the numerator that com- 
pensates the zero in the denominator at k = zk„- Correspondingly, 

1 ^ \r 



propSl 



prop9 



proplO 



propll 



propl2 



gA(a;,x',k) = — y $(x,iK„)*„(x')cot-i^^-^, (4.17) I propl3 



n=l 



where we omitted a term proportional to the l.h.s. of (2.31). This relation is explicitly 
discontinuous at every k — inn and only at these points. Precisely, 

\ jjc ^ 

t/A(a;,x',k) = — $(x, iK„)*„(x')cot"^ — ^ — ^- 

iTT I ksR 

M 

y $(x,iK„)*„(x')+o(l), k-m„, n = l,...,7V. (4.18) 



rn— n+l 



propl4 
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Let us also mention that thanks to (4.15)-(4.17) 

g(x,x','k) |k,^^o= g^{x,x') - -^sgn(K„ -kQ)$(a;, ««„)*„ (x'), 

ri=l 

where it is assumed that k^ y^ k„, n — 1, . . . ,7V. Thanks to (2.31) and (4.13) this gives 

g{x,x',k) \^^^^=g-{x,x\k^^). (4.19) 

On the other side, in vicinity of the points iK„ we have 

g{x, x', k) = g+{x, x') + — $(x, iK„)*„(x') cot-1 f """" ^ 

it: I ksR I 



Af 



y^ $(a;, iKm)*„(x') + o(l), k - ««„, n = l,...,A/'. 



(4.20) 



m—n-\-l 



In a forthcoming pubhcation we wiU show that these properties of the Green's functions 
enable the solution of the problem formulated in the Introduction, i.e., to develop the 1ST 
for the case of perturbed multisoliton potentials (1.3) of the heat operator (1.2). 
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